New cross-phase modulated localized solitons in coupled atomic-molecular BEC 
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The interacting atom-molecule BEC (AMBEC) dynamics is investigated in the mean field ap- 
proach. The presence of atom-atom, atom-molecule and molecule-molecule interactions, coupled 
with a characteristically different interaction representing atom-molecule interconversion, endows 
this system with nonlinearities, which differ significantly from the standard Gross-Pitaevskii (GP) 
equation. Exact localized solutions are found to belong to two distinct classes. The first ones are 
analogous to the soliton solutions of the weakly coupled GP equation, whereas the second non- 
equivalent class is related to the solitons of the strongly coupled BEC. Distinct parameter domains 
characterize these solitons, some of which are analogous to the complex profile Bloch solitons in 
magnetic systems. These localized solutions are found to represent a variety of phenomena, which 
include co-existence of both atom-molecule complex and miscible-immiscible phases. Numerical sta- 
bility is explicitly checked, as also the stability analysis based on the study of quantum fiuctuations 
around our solutions. We also find out the domain of modulation instability in this system. 

PACS numbers: 03.75.Lm , 03.75.Kk 



I. INTRODUCTION 



Molecular BECs have been experimentally realized in 
recent times [iHl]. Co-existence and inter conversion of 
atomic and molecular Bose-Einstein condensates (BECs) 
have been observed experimentally. Raman photoasso- 
ciation is an important process by which the molecular 
species in an AMBEC can be formed. This was investi- 
gated theoretically in Refs pTMl4| . The mean field de- 
scription of the same involves generalization of the Gross- 
Pitaevskii equation to take into account atom-molecular 
two-body scattering, as well as the atom-molecule inter- 
conversion. There are theoretical predictions that the 
atomic and molecular species can show distinct collec- 
tive oscillations in an AMBEC [THS]. In cigar-shaped 
BEC, this rich dynamical system paves the way for ob- 
servation of novel solitons and nonlinear periodic waves, 
akin to the fundamental dark and bright solitons of the 
atomic BEC, in the repulsive and attractive regimes [TSl - 
[T51 [3T] . The fact that GP equation in one dimension 
is the integrable nonlinear Schrodinger equation, which 
admits soliton solutions, has led to considerable theoret- 
ical and experimental investigations of the cigar-shaped 
BEC. Dark solitons, bright solitons and soliton trains 
have been experimentally observed [T51118j . Novel insta- 
bility mechanisms have been proposed for the break up of 
bright soliton and formation of soliton trains, since mod- 
ulation instability has not been adequate in explaining 
the same [HITn]. The two-component BEC (TBEC) is, 
in general non-integrable, having close connection with 
the integrable Manakov system [TO]. The soliton solu- 
tions and their structure and stability has been exten- 
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sively studied for this system, both analytically and nu- 
merically. The mean field equations describing the atom- 
molecule BEC complex, has close similarity, with both 
weakly and strongly coupled atomic BEC. The two-body 
atom-atom, molecule-molecule and atom-molecule scat- 
tering terms are analogous to cubic nonlinearity of the 
standard GP equation, whereas the quadratic nonlinear 
terms arising from atom-molecule conversion is identi- 
cal to the nonlinear interaction term in strongly-coupled 
BEC, in one-dimension |22j. This dual structure of the 
interaction terms, provides a novel form of cross-phase 
modulation, not possible in the conventional TBEC case. 



II. THE MEAN FIELD DESCRIPTION OF 
AMBEC 

In the absence of the trap, the mean-field dynamics 
of the cigar-shaped AMBEC complex is governed by the 
mean-field equations ~~ 
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where N is the total number of atoms in the system. 
Here, e is the binding energy and the terms with coeffi- 
cients Aa, Am and \am dcnotc the effect of atom-atom, 
molecule-molecule and atom-molecule collisions, respec- 
tively. The interaction involving a denote the conversion 
of atoms to molecule and vice-versa. In comparison, the 
mean-field GP equations in the weak and strong coupling 
sectors are given by P^ . 
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under the condition 2aN\<j)\'^ ^ 1 and 



This ansatz solution represents the physical scenario, 
where the asymptotic condensate density of atoms van- 
ishes, and that of molecules reaches a constant value. 
The equations of motion yield the following consistency 
relations: 
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N is 



is the trapping frequency in the radial direction 
the number of atoms in the condensate, a is the scat- 
tering length and ctq is the equilibrium density of atoms, 
far away from the axis. These two arise from the non- 



polynomial interaction term 



when a dimen- 



y/l+2aN\4,\ 

sional reduction of the 3D GP equation is carried out for 
the cigar shaped BEG. In the weak coupling case, one can 
neglect the term in the denominator, yielding the fa- 
miliar cubic nonlinearity, whereas for the strong coupling 

case one gets ^-^4> — \4'\4>- The structure of the solu- 
tions are quite different in these two sectors. The dark 
and bright solitons are of the type A tanh {x — ut) 



and B sech 



(x — ut) 



cos t 



whereas in the later case, 
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the type A sech 



(x — ut) 



for the weakly coupled case, 
one finds the solutions are of 
I- constant [Ml. 



In the following, we highlight the above mentioned sim- 
ilarities between atom-molecular BEG and weak-strong 
coupled atomic BEG. The different nature of the local- 
ized modes in the two different regimes, is pointed out. 
Subsequently, we exhibit the exact solutions of this dy- 
namical system, which realize a new cross-phase mod- 
ulation, arising due to atom-molecule interaction. We 
check the stability of our solution numerically as well as 
based on the study of quantum fluctuations around our 
solutions. We also find out the domain of modulation 
instability. We then conclude with directions for future 
investigation in this rich dynamical system. 



III. SOLITON SOLUTIONS 
A. Sech- Tanh complex soliton pair 

For identifying exact solitonic solutions to the equa- 
tions [1] and [2], we start with the following ansatz for 
the mean fields: 



(j)a{x,t) = acosflsech 
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exp[i (px — fit)], 
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A lengthy but straightforward calculation leads to the 
following solutions: 
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where 



71 — (AaAm - A^„) , 

72 = (2Aa + 2\m - 5Aam) 
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FIG. 1; (top) Bright soliton (5) and (bottom) grey soli- 
ton (6) propagation snapshots at different times. Aa = 
3x10*, Am = 10^A™ = 80.5xl0^a = 20, N = 100, e = 
0,a = 0.000132,6 = 0.000138,6) = 0.4336, ^ = 3.147917, f7 = 
0.094963 and u = 0.122704. 
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FIG. 2: Depiction of for the second class of solution for 
/3 = 1 (top) and for /3 = 3 (bottom). Parameters used are p = 
1, e — 100, N — 100 and a = 41. \4>m\^ has the same profile. 
The figure represents snapshots at discrete times. 



One can now put any valid value of momentum p to 
the above equations and obtain the remaining variables 
exactly in terms of the parameters of the equations (1) 
and (2). From the above equations, it is easy to see, for 
p = 0, that the velocity of the solitons and the depth of 
the grey soliton, are completely fixed by the parameters 
in the mean field equations (1) and (2). Hence, in order 
to obtain grey solitons of varying depths, one needs to 
tune the parameters in the theory. 



The consistency conditions allow only two discrete val- 
ues of the parameter a: 
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B. New class of Soliton solutions 



and 



Due to the nature of the interconversion term in equa- 
tions (1) and (2), only very special classes of solutions 
are allowed in the system. Another class of solutions for 
the AMBEC is presented, which resemble the solutions 
in the strongly coupled BECs [24] . This is a special class 
of solution since the atomic and molecular phases have 
the same density profiles at all times and therefore the 
atoms and the molecules are always in a miscible phase. 
The soliton profiles are given by, 

(f>aix,t) = a [l — cr tanh'^ [7(2:: — ut)]] exp[i {px — ilt)] , 

(19) 

and 

(t>m{x,t) = a [1 — crtanh^ [7(2; — ut)\\ exp[2i [px — Vlt)] 

(20) 



IV. 
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STABILITY UNDER QUANTUM 
FLUCTUATIONS 



(25) 



We now investigate the stability of the obtained solu- 
tions, using the method of C. K. Law et. al. This 
analysis revealed the regime of coupling parameters in 
the theory, in which, the ground state solutions to the 
coupled NLSE, were stable under vacuum fluctuations. 
The above authors identified an eigenvalue, associated 
with the system as the determiner of condensate stabil- 
ity. It plays the same role as the sign of scattering length 
in a single species condensate. We will carry out a similar 
analysis here to find out the parameter domains, where 
the complex bright-grey pair of soliton solution is stable. 
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To perform the analysis, one starts with the second quan- 
tized grand canonical Hamiltonian of the atom-molecule 
BEC: 



At the temperature T = if, we linearize this Hamil- 
tonian by assuming 



tions. The bright-grey pair solution, (5) and (6), that was 
obtained for AMBEC and was tested for its stability. The 
fluctuation part of i^aim) is described by i/'a(m) which 
obey the usual equal-time commutation relations 



[^.(rl, V,(r-')] = 0, = S,Ar~ ^ • (28) 

One obtains the linearized Hamiltonian, by discarding 
terms beyond the second order in the fluctuations. The 
contribution of the fluctuation part to the above Hamil- 
tonian is given by. 



Ipm = 4>ni + i^n 
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where (/)a and (pm are solutions to the coupled GP equa- 
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Taking cue from .25], the AMBEC system is stable if 
all the eigenvalues of M are non- negative, i.e., if M is 
semi- positive. The system is unstable if the lowest eigen- 
value of M is negative. This is justifled by the fact that 
arbitrary fluctuations increasing the energy of the sys- 
tem, should mean that the system was stable in the first 
place. Similarly, if the fluctuations decrease the energy 
of the system, the system was not stable to start with. 
The analysis assumes nothing regarding the nature of the 
fluctuations. The lowest eigenvalue determines the sta- 
bility of the mean fields. It has been pointed out that the 
fluctuations can also be in particle numbers and that the 
consideration of these fluctuations would provide a way 
to probe the effective interactions between particles. 

Evolution of the mean fields is governed by the NLSE 
and small perturbations of the mean fields will remain 
bounded if all the normal mode frequencies (or collective 
excitation frequencies) of the linearized system are real. 
The collective excitation frequencies are defined by 
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where, 



(31) 



with (uifc, W2fc, fifc, f2fc) being the mode functions. A 
semipositive M does indeed guarantee stability of the 
mean fields. A necessary but not sufficient condition for 
mean fields to be unstable (have complex normal mode 
frequencies), is that the lowest eigenvalue of M be nega- 
tive [27]. The final task is to compute the lowest eigen- 
value of M, for varying values of equation parameters. It 
is done by discretizing in space. We fixed time to be 
and assumed that the eigenfunctions of the matrix van- 
ish at ±cxi. We varied Xam, the coupling constant for 
the interaction between atomic and molecular species, 
the particle number N and the interconversion coupling 
constant a. 

From this analysis, we have identified a region in the 
parameter space where the solution makes a transition 
from 'unstable' to 'stable'. Fig (3) shows the least eigen- 
value of M, computed as a function of atom-moleculc in- 
teraction coefficient Xam- Different curves were obtained 
from using different values of the interconversion coeffi- 
cient a. We then use the stable parameter domain in 
simulating our solution. 
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FIG. 3: Stable parameter domain under quantum fluctua- 
tions to the coupled soliton solutions of the bright-grey type. 
Xa = SxlO"*, A„ = 10^ e = and N=100. Varying the par- 
ticle number (N) does not change the plots by a considerable 
extent. The curves are very sensitive to the interconversion 
coefRcient a. The 'stable-unstable' transition occurs at dif- 
ferent values of Xam for different a. 



V. NUMERICAL STUDY OF SOLITON 
SOLUTIONS 

We have further studied the evolution of the bright- 
grey soliton solution numerically. The following are the 
results of the simulations carried out on the bright-grey 
soliton pair. Figs 4 and 5 were obtained by the stable 
finite difference scheme, Crank Nicolson (CN). It is an 
implicit scheme and accumulates second order error in 
both space and time steps. 

The following simulations were implemented using 
the coupled split-step and CN scheme 28j. The so- 
lution is first evolved using split-step method, using 
only the terms containing atom-atom, molecule-molecule 
and atom-molecule interactions and the binding energy. 
The CN algorithm next evolves this evolved part, using 
the atom-molecule interconversion and dispersion terms. 
The results are shown in the following figures. 

The results depicted in Figs 4 and 5 suggest that the 
solutions are stable when the parameters are chosen from 
the stable region of Fig 3 and unstable when the least 
eigenvalue is highly negative. However, an exhaustive 
study needs to be done to conclude if the analysis per- 
formed on our solution correctly predicts stable- unstable 
domains. It would also be worthwhile to see what predic- 
tions linear stability analysis will have for stable- unstable 
domains and whether or not there is any overlap in the 
domains predicted by these two analyses. 




FIG. 4: Results of numerical evolution of the bright and grey 
pair solitons in the stable domain. The parameters used were 
obtained from the stable domain of Xam in Fig (3), and are 
A,, = 30xl0^ A™ = 10^ Xam = 80.5xl0^ a = 20, N= 100, 
and e = 0. The least eigenvalue in this case is -1-0.0466, indi- 
cating stable domain. 

are generated, is through modulation instability (MI). In 
this case, the continuous wave solution becomes unstable. 
MI of a nonuniform initial state in the presence of a har- 
monic potential has been studied both analytically and 
numerically in the context of the mean field of the BEC 
[30] . The analysis of MI in AMBEC is similar to that in 
the two component BEC (TBEC), but not exactly the 
same [?T] . 

A. Gain Spectrum 

To find out the domain of MI in any system, in gen- 
eral one proceeds as follows. We first find out continuous 
wave solutions to the coupled GP equations that are fixed 
in space. Subsequently one applies space-time dependent 
perturbation to this solution and finally, the gain spec- 
trum is obtained. For this purpose, we use the ansatz 



VI. MODULATION INSTABILITY 



(fiaix, t) = (0qo + ea{x, t)) exp{iipat) , (32) 

and 



We now proceed to study the possibility of modula- 
tional instability in this system. In NLSE, the standard 

way in which bright solitons and solitary wave structures (t'm{x, t) = {(l>mO + em(a;, t)) ex.p{iipmt) ■ (33) 
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Then we assume, 




ea{x, t) = Ua exp[i(fca; — ^It)] + w* cxp[— i(fc.T — ^It)] 

(34) 

Cmix, t) ~ Ujn exp[i(fcx — Vlt)] + v^^ exp[— i(fca; — Vlt)] 

(35) 



FIG. 5: Results of numerical evolution of the bright and grey 
pair solitons in the unstable domain. The parameters used 
are A = 0.15, A™ = 0.50, A^™ = 0.506, a = 41.5, N=1000, 
and e = —900. The least eigenvalue in this case is —4475.8 
indicating unstable domain. 



where Ua(m) i'*^^-) are to be determined. The fol- 
lowing represents the consistency condition, that (32)- 
(35) yield valid solutions to the coupled GP equations, 
in terms of a matrix determinant: 



71 (fc) 





(36) 



where 

71 = e/2 + Z\aN4>l^ + \amN4>l,a + aV^N (p^^ (37) 

72 = + \aN4>l^ + XarnN(g,^ - aV2N (f^mO (38) 

73 = A:V4 + e + 3A™iV0^o + AamA^<^aO (39) 

One obtains a quadratic equation in fi^, which gives 
two roots. MI sets in when fl^ < 0. The growth rate is 
given by the imaginary part of il. Plotting this as a func- 
tion of varying parameters, one gets the gain spectrum. 
We study one of the branches (fi^ branch), correspond- 
ing to one of the roots of the above mentioned quadratic 
equation. The other branch U~ was also studied and 
yielded a similar gain spectrum. 

Fig 7 gives the gain spectrum for MI in the AMBEC 
system. The two modulationally unstable species, atoms 



and molecules, may appear as propagating periodic or 
localized solitary waves. 



VII. CONCLUSION 

In conclusion, we have found new cross-phase modu- 
lated localized soliton solutions for an AMBEC. Owing 
to the difference in nature of the atom-molecule inter- 
conversion terms in Eqs (1) and (2) and the equivalent 
terms in a TBEC, the solutions vary significantly in these 
two cases. Many of the solutions valid for the case of a 
TBEC, do not obey Eqs (1) and (2). We have iden- 
tified three solutions for the case of AMBEC and have 
devised a mechanism for obtaining more. We have an- 
alyzed one of these solutions, the bright-grey pair, for 
stability under quantum fluctuations, by performing an 
analysis presented in 125] . The analysis helped us predict 





FIG. 7: Gain Spectrum for MI in AMBEG. The following are 
the values of the parameters that were used: = 30x10'^, 
\^ = 10^ Aarn = 90xl0^ a = 0, N = 100, e = 0. 



FIG. 6: Results of numerical evolution of the bright and grey 
pair solitons in the stable domain. The parameters used were 
obtained from Fig 3 and are Aa = 30x10^, Am ~ 10^, A^m = 
80.5x10^, a = 20, N= 100, and e = 0. Red and blue denote 
the numerically evolved solutions, black and green, denote 
analytical solutions (5) and (6). A good match between the 
two indicates that the solutions in this domain are stable. 



stable and unstable regions in the parameter space. This 
was supported by the numerical simulations (Figs. 4-6). 
We also obtained the domain of modulation instability 
in the AMBEC system. 
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